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摘 要 : 研究 山地 过 渡 带 植被 物候 格局 以 及 植被 物候 变化 的 驱动 机 制 具有 重要 意 


点 实验 室 ,陕西 西安 710127) 


基于 2000 一 2018 年 MODIS 


NDVI 遥 感 数据 ,利用 双 Logistic 曲线 拟 合法 提取 2000 一 2018 年 秦岭 山地 植被 物候 参数 ,结合 气温 和 降水 数据 ,运用 
趋势 分 析 、 相 关 分 析 .元 余 分 析 等 方法 ,研究 秦岭 山地 不 同 物候 期 的 时 空 变化 特征 和 对 气象 因子 的 响应 。 结 果 表 
明 :(1) 秦岭 山地 物候 始 期 集中 于 60~102 d ,物候 未 期 主要 集中 在 315~345 d, 生 长 季 长 度 集中 在 225~255 d ,空间 上 


看 ,具有 明显 的 垂直 地 带 分 布 特征 , 随 海拔 升 高 ,物候 始 


期 ,未 期 .生长 期 长 度 分 别 呈 延迟 ,提前 和 缩短 趋势 (2) 从 


植被 物候 的 年 际 变化 来 看 ,物候 始 期 62.25% 的 区 域 呈 现 提前 趋势 ,53.42% 的 区 域 物 候 末期 呈现 推迟 趋势 ,其 中 


59.18% 的 区 域 植被 生长 期 长 度 延 长 , 仅 有 5.12% 位 于 中 东部 中 高 海拔 区 的 部 分 区 域 显 著 延 长 。(3) 植被 物候 始 期 大 
部 分 区 域 与 2 一 $ 月 平均 气温 和 降水 呈 负 相关 关系 ,分 别 占 总 面积 的 50.29% 和 65.24% ;物候 末期 与 8 一 11 月 平均 气 


温和 月 降水 量 主要 呈正 相关 关系 ,分 别 占 总 区 域 的 66.63% 和 59.77%;(4) 宛 余 分 析 (Redundancy analysis ,RDA) 结 果 


表明 秦岭 山地 植被 物候 受 当 季 和 六 
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植被 物候 是 植物 长 期 适应 环境 的 季节 性 变化 
而 形成 的 生长 发 育 节律 一, 植物 生长 物候 不 仅 能 
反映 季节 变化 ,还 对 生态 系统 碳水 以 及 能 量 的 循 
环 具有 决定 性 作用 中。 研究 植被 物候 的 变化 规律 
及 其 对 气象 因子 的 响应 程度 ,有 助 于 理解 和 完善 植 
被 一 气候 关系 。 

基于 高 时 空 分辨 率 遥感 数据 反 演 的 NDVIEVI 
等 植被 指数 ,国内 外 学 者 对 全 球 不 同 尺 度 物 候 监 测 
进行 了 研究 ,探讨 了 不 同 区 域 植被 物候 的 变化 及 其 
对 气象 因子 的 响应 。 欧 亚 大 陆 北部 1982—2004 4 
遥感 数据 分 析 植被 物候 变化 结果 显示 生长 期 开始 日 
期 平均 提前 了 8 d, 结 束 日 期 则 平均 推迟 了 3.6 d“; 
在 我 国 的 东北 5 西北 青藏 高 原 "” 等 地 的 植被 物 
候 也 呈现 出 春季 物候 提前 ,秋季 物候 期 推迟 ,生长 
季 延 长 的 现象 ,并 且 气 温 在 物候 变化 中 起 主导 作 
用 ;此 外 ,降水 也 是 影响 植被 物候 变化 不 可 缺少 的 


1 期 气象 因子 的 共同 影响 ,春季 物候 受 物候 发 生 期 的 影响 大 于 上 年 冬季 ,物候 来 
期 前 期 的 气温 和 降水 对 秋季 物候 的 影响 强度 大 于 秋季 ,不 同 坡 向 物候 变化 对 气象 因子 的 响应 程度 有 所 差异 。 


区 遥感 植被 指数 峰值 出 现在 干旱 季节 后 期 至 雨季 
前 期 ; 钥 佳 星 等 % 利用 CIMMS AVHRR 遥感 数据 分 
析 我 国 东北 地 区 物候 对 气候 的 响应 ,发 现 植 被 秋季 
物候 变化 主要 受降 水 影响 ;以 上 的 研究 可 以 看 出 ， 
植被 物候 变化 存在 区 域 性 差异 ,同一 地 区 时 间 范 
围 、 物 候 期 提取 所 选 方法 和 遥感 数据 源 不 同 ,植物 
物候 表现 也 有 所 差异 。 因 此 ,针对 特定 的 研究 区 
域 , 选 择 合适 的 时 间 序 列 和 物候 期 提取 方法 ,精确 
人 研究 物候 变化 对 气象 因子 的 响应 至 关 重要 。 

秦岭 一 大 巴山 是 我 国 南北 过 渡 带 的 主体 ,是 控 
制 中国 地 理 格局 和 空间 变化 的 枢纽 ,秦岭 作为 我 
国 重 要 的 地 理 生 物 分 界线 六 ,垂直 地 带 性 显著 A 
然 地 理 跨度 大 ,植被 类 型 复杂 多 样 ,生态 系统 较为 脆 
弱 。 同 时 随 着 全 球 气 候 变 化 ,秦岭 气温 降水 ”等 
气象 因子 均 产 生 了 显著 变化 。 随 着 气候 变化 秦岭 
山地 的 植被 物候 会 产生 怎样 的 变化 ” 物候 期 趋势 
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关系 如 何 ?为 了 解决 这 些 问 题 , 本 文 利用 双 Logistic 
拟 合 ,结合 NDVI 数 据 和 地 面 气温 、 降 水 资料 ,提取 
T 2000—2018 年 秦岭 山地 植被 物候 生长 季 始 期 ,未 
期 和 生长 期 长 度 , 分 析 秦 岭 山地 植被 物候 变化 趋势 
和 规律 ,探讨 物候 时 空 变化 及 其 对 气象 因子 的 响应 
程度 ,识别 影响 秦岭 山地 物候 变化 的 主要 范围 时 
段 , 旨 在 为 全 球 变化 背景 下 山地 生态 系统 的 监测 和 
保护 提供 参考 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

人 研究 区 为 陕西 境内 的 秦岭 , 即 “ 陕 西 秦岭 ”, 北 
临 渭河 ,南面 汉江 ,地 理 位 置 为 105°30'~110°05'E， 
32?40'~34?35'N ,包含 宝鸡 西安 .商洛 等 七 个 地 级 
市 ,海拔 高 度 177~3746 m( 图 1)。1959 一 2018 年 平 
均 气 温 12.1~14.4% ,年 降水 量 495~1137 mm。 秦 岭 
是 横 吾 于 中 国 中 部 的 东西 走向 的 巨大 山脉 ,特殊 的 
地 理 位置 和 复杂 的 气候 条 件 使 得 南北 坡 形成 了 差 
异 明 显 的 垂直 带 谱 , 拥 有 常 绿 阔 叶 林 、 针 冰 叶 混交 
林 、 高 山 草 多 等 种 类 繁多 的 植被 类 型 和 褐 土 . 棕 壤 、 
黄 棕 壤 等 多 样 化 的 土壤 类 型 。 
1.2 数据 来 源 与 处 理 
1.2.1 NDVI 数 据 采用 2000 一 2018 年 的 MOD13Q1 
数据 集 ,数据 集 来 源 于 美国 国家 航空 航天 局 NASA 
(https://www.nasa.gov/) ,时 间 分 辩 率 16 d, 空间 分 辩 


106°E 107°E 108°E 


率 为 250 mx250 m。 首 先 使 用 MRT (MODIS Repro- 
jection Tool, MRT) 工 具 从 MOD13Q1 数据 集 提 取 第 
一 波段 NDVI, 进 行 投影 转换 和 图 像 拼 接 , 其 次 采用 
MVC 最 大 合成 法 消除 残留 的 云 .大 气 和 太阳 高 度 
等 ,最 后 利用 研究 区 边界 对 数据 进行 掩 模 提 取 , 获 
得 2000 一 2018 年 的 月 NDVI 影 像 。 

1.2.2 气象 数据 气象 数据 为 研究 区 内 32 个 国家 
标准 气象 台 站 2000 一 2018 年 的 气温 .降水 数据 。 数 
据 来 源 于 中 国 气象 数据 科学 共享 服务 平台 (www. 
geodata.cn) 和 陕西 省 气象 局 。 以 站 点 数据 为 基础 ， 
采用 空间 插值 法 将 点 状 数据 转换 为 面 状 气象 数据 ， 
最 后 利用 研究 区 边界 进行 掩 模 提 取 , 得 到 秦岭 山地 
气温 、 降 水 插值 数据 。 其 中 ,山地 气温 条 件 复 杂 , 采 
用 前 人 人 研究 所 得 的 秦岭 气温 垂直 递减 率 结果 进行 
空间 插值 ,降水 采用 普通 Kriging 插 值 法 。 

1.2.3 其 他 数据 ”DEM 数据 采用 V2 版 本 的 AS- 
TER GDEM 产品 (空间 分 辨 紊 为 30 m) ,来 源 于 地 理 
空间 数据 云 (http://www.gscloud.cn) ,对 高 程 数据 进 
行 投影 转换 、 重 采样 ,最 后 掩 模 提取 ,得 到 空间 分 辨 
率 250 m 的 秦岭 地 形 数据 。 

1.3 研究 方法 

1.3.1 植被 物候 各 感 提 取 方 法 利用 TIMESAT 软 
件 对 植被 物候 期 进行 提取 。 首 先 ,虽然 MOD13Q1 
数据 集 经 过 最 大 值 合 成 法 (MVC) 处 理 降 低 了 云 和 
大 气 的 影响 ,但 太阳 高 度 角 观测 角度 .臭氧 ATA 
胶 等 一 些 随机 性 因素 的 存在 仍 使 得 曲线 季节 变化 
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图 1 研究 区 示意 图 


Fig. 1 Sketch map of study area 
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趋势 不 明显 ,因此 ,本 文选 用 软件 包 中 的 双 Logis- 
tic 拟 合法 (公式 1) 对 2000 一 2018 年 NDVI 影 像 进行 
拟 合 处 理 。 其 次 ,采用 DL 模型 根据 研究 区 提取 
2000—2018 年 的 植被 物候 参数 , 即 生长 季 始 期 
(Start of the Growing Season ,SOS) 生长 季 末 期 (End 
of the Growing Season,EOS) 和 生长 季 长 度 (Length of 
the Growing Season ,LOS)。 
1 1 


t; QQ = = 
alts a J ou 一 a,-t 
1+exp 1+exp = 


Qa, 4 


(1) 
式 中 ;1 表示 时 间 尺 度 ; aa, 表示 形状 和 位 置 参 
数 ; a 为 生长 季 始 期 所 处 位 置 , a, 为 生长 季 始 期 变 
化 率 , a 为 生长 期 末期 所 处 位 置 , a 为 生长 期 末期 
变化 率 。 
1.3.2 趋势 分 析 采用 一 元 线性 回归 分 析 法 逐 像 元 
模拟 2000 一 2018 年 生长 季 始 期 ,未 期 .生长 期 长 度 
的 趋势 斜率 ,公式 如 下 : 
nx YlixP)- Yi 


on (2) 


Slope = 


IWP: n 为 年 数 ; P, 为 第 i 年 的 生长 参数 ; Slope 为 
一 元 线性 回归 方程 的 斜率 , Slope >0 说 明生 长 参数 
总 体 变 化 是 推迟 (延长 ) 的 ,反之 是 提前 (缩短 ) 的 。 
并 对 变化 趋势 进行 显著 7 检验 。 

1.3.3 相关 性 分 析 采用 偏 相 关 分 析 揭 示 气 温和 降 
水 与 植被 物候 变化 之 间 的 响应 关系 ,公式 如 下 : 


(3) 


式 中 : rs 表示 x 与 y 之 间 的 相关 系数 ,表示 两 要 素 
之 间 的 相关 程度 ,+，>0, 表 示 正 相关 ,<0, 表 示 负 
相关 ; Xx、 了 分 别 表示 要 素 样 本 值 的 平均 值 。 

1.3.4 TAAS TCRA HN (Redundancy analysis, 
RDA) 是 多 变量 环境 梯度 分 析 ” ,是 响应 变量 矩阵 与 
解释 变量 之 间 多 元 多 重 线性 回归 的 拟 合 值 矩 阵 分 
析 , 目 前 已 广泛 应 用 于 生物 、 环 境 、 土 壤 等 研究 领 
域 '“*。 本 文选 取 秦 岭 山 地 物候 开始 和 结束 时 间 前 
8 个 月 的 气温 和 降水 作为 气象 因子 ,采用 元 余 分 析 
量化 物候 期 参数 与 气象 因子 之 间 的 关系 。 使 用 
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Canoco 5 软件 进行 元 余 分 析 Be FF KY BE (Monte 
Carlo Test) 和 RDA 排序 图 绘制 ,物种 和 环境 变量 图 
中 ,箭头 所 处 的 象限 表示 环境 因子 与 排序 轴 间 的 正 
负 相 关 性 ,箭头 连 线 越 长 ,表示 某 个 环境 因子 与 物 
种 分 布 的 相关 程度 越 大 ,通过 物种 箭头 垂直 投影 到 
环境 变量 箭头 的 线 上 来 近似 估算 种 类 与 环境 变量 
之 间 的 相关 性 ,投影 点 沿 箭 头 指示 的 方向 越 远 , 相 
关 性 越 高 ,坐标 原点 附近 的 投影 点 表示 相关 性 接近 
E ,投影 点 位 于 反方 向 ,表明 为 负 相关 。 


2 结果 与 分 析 


2.1 平均 物候 期 空间 分 布 

从 2000—2018 年 秦岭 山地 植被 物候 期 (SOS、 
EOS、LOS) 的 空间 分 布 可 以 看 出 (图 2) ,物候 始 期 均 
值 主要 集中 分 布 在 第 60~102 d(3 月 至 4 月 中 旬 ), 占 
研究 区 总 面积 的 81% ,物候 始 期 呈现 从 南北 两 侧 低 
海拔 地 区 向 中 心 高 海拔 区 域 推迟 的 规律 , 北 坡 的 地 
势 低 平 区 , 始 期 早 于 第 60 d, 而 在 一 些 高 海拔 和 山 崔 
线 附 近 的 地 区 , 始 期 则 晚 于 第 114 d; 由 图 2b 可 以 看 
出 ,未 期 北 坡 平均 发 生 于 第 304 d, 南 坡 平均 发 生 于 
第 325 d, 北 坡 较 南 坡 平 均 早 21 d, 其 中 汉中 盆地 和 
安康 月 河川 道 两 地 ,物候 未 期 较 早 ,可 能 是 由 于 两 
地 地 势 低 平 ,农业 比较 发 达 , 受 人 类 生产 生活 影响 
大 ,农作物 生长 期 短 于 和 森林 植被 ,成 熟 期 较 早 ;物候 
生长 期 长 度 , 除 南北 两 侧 少 部 分 低 平 区 ,基本 呈现 
出 从 南北 两 侧 至 山脊 线 逐 渐 缩 短 的 规律 。 总 体 来 
看 ,物候 期 分 布 格局 与 海拔 高 度 和 地 形 保 持 了 较 好 
的 一 致 性 ,海拔 越 高 ,物候 始 期 越 晚 ,未 期 越 早 , 生 
长 期 越 短 。 
2.2 变化 趋势 分 析 

从 秦岭 山地 物候 期 参数 的 年 际 变 化 趋势 可 知 
(图 3)。 物 候 始 期 62.25% 的 区 域 呈 现 提 前 趋势 ， 
37.75% 呈 推迟 趋势 ,但 是 提前 和 推迟 的 程度 较 小 ， 
其 中 显著 提前 的 区 域 仅 占 2.12%(P<0.05) ,推迟 趋 
势 达 到 显著 性 (P<0.1) 的 区 域 仅 占 3.26% ;显著 推迟 
的 区 域 主要 分 布 在 南北 葛 低 海拔 区 和 南部 几 个 受 
人 类 影响 较 大 的 区 域 。 

物候 末期 主要 呈现 推迟 的 趋势 (53.42% ) ,其 
中 ,7 检验 结果 表明 ,大 部 分 区 域 呈 不 显著 推迟 和 不 
显著 提前 趋势 (83.37% ) , 仅 有 7.62% 的 区 域 显著 推 
述 (P<0.05), 主 要 分 布 在 秦岭 北 芳 低 海拔 地 区 ,显著 
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图 2 2000 一 2018 年 秦岭 山地 森林 植被 物候 期 多 年 均值 空间 分 布 


Fig. 2 Spatial distribution of multi-annual mean forest phenology in Qinling Mountains from 2000 to 2018 
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Fig. 3 Spatial distribution and significance of inter-annual changes in phenology of 
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forest vegetation in Qinling Mountains from 2000 to 2018 
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提前 的 区 域 仅 有 2.95% (P<0.05) , 主要 分 布 于 秦岭 。 部 地 区 。 
西 侧 。 物候 未 期 与 8 一 11 月 平均 气温 和 月 降水 量 相关 
秦岭 山地 植被 生长 期 59.18% 的 区 域 呈 延长 趋 ”分析 空间 分 布 可 以 看 出 (图 4c, 图 4d) ,秦岭 山地 植 
势 ,40.82% 的 区 域 呈 缩短 趋势 。 其 中 ,7 检验 结果 显 ”被 物候 未 期 大 部 分 区 域 与 同期 气温 和 降水 呈正 相 
示 , 和 森林 植被 生长 季 长 度 的 变化 不 显著 ,50.15% 和  ” 关 关 系 , 分 别 占 总 面积 的 66.63% 和 59.77%, 显 著 相 
37.46% 的 区 域 分 别 旦 现 不 显著 延长 和 缩短 的 趋势 ” 关 的 面积 分 别 为 5.64% 和 8.27% ,这 表明 温度 升 高 ， 
(P<0.05) ,对 于 变化 显著 的 区 域 而 言 ,1.56% 显 著 缩 。 降水 量 增高 使 得 秋季 物候 期 推迟 ,空间 上 ,表现 为 
短 区 域 (P<0.05 ) 主 要 分 布 在 易 受 人 类 活动 干扰 的 低 。 随 气温 升 高 而 物候 期 推迟 的 区 域 主要 位 于 高 海拔 
海拔 区 ,5.12% 显 著 延 长 域 (P<0.05) 主 要 和 集 在 中 东 ”地 区 ,表明 物候 变化 与 地 形 起 伏 具 有 一 致 性 。 


部 中 高 海拔 区 。 2.4 TAD 
2.3 物候 与 气象 因子 相关 性 分 析 由 于 气候 变化 和 物候 变化 存在 时 间 上 的 滞后 


从 物候 始 期 与 2 一 5 月 平均 气温 和 降水 的 空间 ”性 外 , 且 已 有 研究 表明 物候 期 不 仅 受 当 季 的 气温 影 
相关 分 析 可 以 看 出 (图 4a, 图 4b) ,秦岭 山地 植被 物 响 ,同时 前 期 气温 与 物候 期 具有 显著 相关 关系 "。 
修 始 期 大 部 分 区 域 与 同期 气温 和 降水 呈 负 相关 关 ”因此 ,本 文选 取 物 候 始 期 和 末期 前 8 个 月 的 气温 、 降 
系 , 分 别 占 总 面积 的 50.29% 和 65.24% ,从 显著 性 水 ”水 作为 气象 因子 ,应 用 宛 余 分 析 (RDA) 量 化 秦岭 山 
平 来 看 ,植被 物候 始 期 与 气温 显著 的 区 域 最 多 地 南北 两 侧 物候 期 参数 与 气象 因子 之 间 的 关系 。 
(12.43% ) ,这 表明 气温 降低 .降水 减少 会 使 得 植被 ”图 5$~ 图 7 中 ,向 量 越 长 表示 气象 因子 越 重要 ,气象 向 
生长 季 始 期 推 后 。 空 间 上 ,2 一 5 月 平均 气温 对 西部 ” 量 与 物候 向 量 夹 角 的 余弦 值 代表 之 间 的 相关 系 
和 北部 大 部 分 区 域 的 物候 始 期 具有 负面 影响 ,平均 ” 数 。 二 者 方向 相同 表明 具有 较 强 正 相 关 关系 ,方向 
降水 和 物候 始 期 在 北 侧 低 海 拔 区 多 为 正 相 关 关 系 ， ”相反 表明 具有 较 强 负 相 关 关 系 , 垂 直 表 示 不 相关 。 
其 余地 区 相关 系数 正 负 值 分 布 均匀 ,这 可 能 是 因为 由 图 5 可 以 看 出 ,所 选 的 气象 因子 可 以 解释 
西部 地 区 海拔 较 高 ,而 北部 地 区 纬度 偏 北 ,2 一 5 月 SOS 和 EOS 的 变化 特征 ,不同 气象 因子 前 两 轴 累 积 
气温 回升 速度 较 慢 ,使 得 物候 始 期 时 间 晚 于 秦岭 南 。 解释 的 信息 量 都 为 1, 表 明 排 序 轴 反映 了 气象 因子 
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图 4 秦岭 山地 植被 物候 期 对 气温 和 降水 响应 的 空间 分 布 


Fig.4 Spatial distribution of response of vegetation phenological period to temperature and precipitation in Qinling Mountains 
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(a) 物候 始 期 (SOS) (b) 物候 未 期 (EOS) 
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注 :数字 表示 月 份 ,7 表示 气温 ,P 表 示 降 水 
图 5 秦岭 植被 物候 期 与 气象 因子 的 RDA 排 序 网 
Fig. 5 RDA ordination diagram of the vegetation phenology 


period and meteorological factors in Qinling Mountains 
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注 :数字 表示 月 份 ,P 为 前 一 年 数据 ,C 为 当年 数据 。 
图 6 秦岭 物候 始 期 与 气象 因子 的 RDA 排 序 图 
Fig.6 RDA ordination diagram of the SOS and 


meteorological factors in Qinling 
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图 7 秦岭 物候 末期 与 气象 因子 的 RDA 排 序 图 
Fig.7 RDA ordination diagram of the EOS and 


meteorological factors in Qinling 


与 物候 期 的 绝 大 部 分 信息 。 春 季 物 候 发 生 期 2 一 5 
月 的 气温 对 物候 始 期 的 影响 最 强 ,与 南 坡 和 北 坡 物 
候 始 期 均 为 负 相关 关系 ;2 一 5 月 的 降水 与 北 坡 的 物 


候 始 期 呈正 相关 关系 ,与 南 坡 物 候 始 期 相关 性 不 显 
著 ; 上 年 10 月 至 当年 1 月 对 物候 始 期 的 影响 最 小 。 
秋季 物候 发 生前 期 4 一 7 月 降水 对 物候 末期 的 影响 
最 大 ,与 南北 坡 物候 末期 均 呈 显著 正 相 关 关 系 ;4 一 
7 月 的 气温 对 物候 末期 的 发 生 影响 次 之 ;8 一 11 月 的 
气温 与 南北 坡 物候 末期 均 呈 正 相 关 关 系 ,8 一 11 月 
降水 与 南 坡 物候 末期 呈正 相关 关系 ,与 北 坡 物候 末 
期 基本 不 相关 。 综 上 所 述 , 秦 岭 山地 物候 发 生 不 仅 
受 物候 期 气象 因子 的 影响 ,物候 发 生 期 前 期 的 气象 
因子 影响 同样 重要 ;春季 物候 受 物候 发 生 期 的 影响 
大 于 上 年 冬季 , 而 物候 末期 前 期 的 气温 和 降水 对 秋 
季 物 候 的 影响 强度 大 于 秋季 。 

为 了 识别 影响 植被 物候 的 具体 时 段 ,针对 气温 
和 降水 要 素 单独 出 图 。 由 于 每 次 分 析 仅 选 用 一 个 
物候 参数 ,使 得 第 一 轴 和 第 二 轴 相 关系 数 为 零 ， 
此 RDA 排序 图 中 仅 可 以 通过 象限 判断 物候 期 与 气 
象 因子 的 正 负 相关 性 ,通过 向 量 长 短 识别 影响 物候 
期 最 显著 的 月 份 。 从 图 6 和 表 1 可 以 看 出 ,秦岭 北 
坡 ,物候 始 期 与 上 年 10 一 12 月 的 平均 气温 相关 性 
低 , 绝 对 值 小 于 0.49 ,表明 前 一 年 气温 因子 对 于 SOS 
影响 不 大 ;物候 始 期 主要 受 当 年 4 月 的 气温 控制 , 相 
关系 数 达 到 0.78; 前 一 年 12 月 的 降水 量 与 SOS 呈 显 
著 正 相关 关系 ,表明 物候 始 期 的 发 生 受 上 年 冬季 的 
降水 量 影响 大 于 当年 春季 的 降水 。 秦 岭南 坡 ,物候 
始 期 与 当年 3 月 .4 月 的 气温 呈 显 著 负 相关 关系 , 表 
明 气 温 越 高 南 坡 春季 物候 期 起 始 时 间 越 早 ; 物候 始 
期 与 上 年 11 月 的 降水 呈 显 著 负 相关 关系 ,表明 冬季 
降水 越 多 不 利于 植物 的 生长 发 芽 。 由 图 6 可 以 看 
出 ,秦岭 北 坡 物候 始 期 主要 受 当年 4 月 气温 和 上 年 
12 月 降水 影响 ,而 秦岭 南 坡 主要 受 当年 3 一 4 月 气温 
和 上 年 11 月 降水 影响 ,这 可 能 与 秦岭 山地 位 于 我 国 
暧 温带 与 亚热带 分 界线 有 关 , 南 坡 向 阳 , 北 坡 背 阴 ， 
二 者 所 受 太阳 辐射 不 同 , 植 被 物候 变化 对 气象 因子 
的 响应 也 有 所 差异 。 

从 图 7 和 表 2 可 以 看 出 ,秦岭 北 坡 , 物 候 末 期 与 
9 月 .10 月 气温 呈 较 明显 的 负 相 关 关 系 , 相 关系 数 
在 -0.53~-0.41 ,与 物候 前 期 的 4 一 6 月 的 气温 呈 较 明 
显 的 正 相 关 关 系 , 其 余 月 份 的 影响 较 小 ;其 中 ,4 月 
降水 与 EOS 旺 显著 正 相关 关系 ,表明 前 期 的 降水 对 
秋季 物候 有 显著 影响 。 秦 岭南 坡 , 物 候 期 与 7 月 气 
温 呈 显著 正 相 关 关 系 , 与 10 月 呈 显 著 负 相关 关系 ; 
南 坡 5 月 的 降水 对 物候 末期 起 显著 的 正面 影响 , 相 
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表 1 2000 一 2018 年 秦岭 山地 物候 始 期 RDA 排 序 前 两 轴 与 气象 因子 的 相关 系数 


Tab.1 Correlation of two SOS's RDA axes and meteorological factors of each plot in Qinling Mountains from 2000 to 2018 


AG 北 坡 气温 北 坡 降水 南 坡 气 温 南 坡 降水 
第 1 轴 排 序 ”第 2 轴 排 序 。” 第 1 轴 排 序 ”第 2 轴 排 序 ” 第 1 轴 排 序 ”第 2 轴 排 序 ” 第 1 轴 排 序 ”第 2 轴 排 序 

10 -0.30 0.00 -0.04 0.00 0.39 0.00 0.50 0.00 
11 0.49 0.00 0.26 0.00 -0.14 0.00 -0.82 0.00 
12 0.29 0.00 0.69" 0.00 0.30 0.00 0.03 0.00 
1 -0.01 0.00 -0.04 0.00 0.20 0.00 0.61 0.00 
2 -0.14 0.00 -0.53 0.00 -0.03 0.00 0.05 0.00 
3 -0.29 0.00 0.01 0.00 -0.71 0.00 -0.19 0.00 
4 0.78" 0.00 -0.56 0.00 -0.49° 0.00 0.34 0.00 
5 0.15 0.00 -0.12 0.00 0.20 0.00 0.11 0.00 


注 :* 表 示 通 过 0.05 可 信 度 检验 的 显著 因子 。 下 同 。 


表 2 2000 一 2018 年 秦岭 山地 物候 末期 RDA 排 序 前 两 轴 与 气象 因子 的 相关 系数 


Tab. 2 Correlation of two EOS's RDA axes and meteorological factors of each plot in Qinling Mountains from 2000 to 2018 


AG 北 坡 气温 北 坡 降水 南 坡 气 温 南 坡 降水 
第 1 轴 排 序 ”第 2 轴 排 序 。” 第 1 轴 排 序 ”第 2 轴 排 序 第 1 轴 排 序 ”第 2 轴 排 序 ”第 1 轴 排 序 ”第 2 轴 排 序 

4 0.52 0.00 0.82° 0.00 0.13 0.00 0.03 0.00 

5 0.47 0.00 0.20 0.00 0.16 0.00 0.80° 0.00 

6 0.41 0.00 0.06 0.00 -0.01 0.00 -0.61 0.00 

7 -0.29 0.00 -0.38 0.00 0.27 0.00 0.15 0.00 

8 -0.26 0.00 -0.20 0.00 0.15 0.00 0.23 0.00 

9 -0.41 0.00 0.10 0.00 -0.27 0.00 -0.04 0.00 

10 -0.53 0.00 -0.18 0.00 -0.64 0.00 -0.27 0.00 

11 0.03 0.00 0.46 0.00 0.01 0.00 0.40 0.00 


关系 数 达 到 0.80。 秦 岭 山 地 南北 两 侧 物 候 未 期 存在 
差异 , 北 侧 植被 生长 结束 时 间 主 要 受 当年 9 一 10 月 
气温 ,4 月 降水 影响 ,而 南 坡 主要 受 当 年 9 一 10 月 气 
温和 5 月 降水 影响 。 


3 讨论 


地 表 植 被 是 区 域 生态 健康 的 鲜明 指示 因子 ,使 
用 宏观 遥感 数据 可 以 实现 对 森林 植被 的 动态 监测 ， 
反 演 山地 生态 系统 植物 物候 的 变化 情况 ,分析 其 对 
气象 因子 的 响应 程度 。 通 过 分 析 ,发现 秦岭 山地 
2000 一 2018 年 植被 物候 主要 呈现 物候 始 期 提前 , 物 
候 末期 推迟 ,生长 期 延长 的 趋势 ,这 与 其 他 学 者 的 
DEA — Be 

研究 结果 表明 ,秦岭 山地 植被 物候 期 受气 温和 
降水 的 共同 作用 ,植被 物候 变化 不 仅 受 物 候 发 生 期 
当 季 的 气象 因子 影响 ,同时 物候 前 期 的 气象 因子 也 
会 产生 作用 ,并 且 秦 岭 山 地 南北 两 侧 物候 变化 对 气 


象 因子 的 响应 存在 差异 。 物 候 始 期 秦岭 北 坡 与 当 
年 2 一 $ 月 的 气温 呈 显 著 正 相关 关系 ,温度 升 高 , 太 
阳 辐 射 增强 ,有 效 积温 提前 来 到 ,从 而 使 得 春季 
物候 期 提前 ;秦岭 南北 两 侧 物 候 始 期 均 受 上 年 10 月 
至 当年 1 月 降水 的 控制 , 呈 负 相关 关系 ,这 可 能 是 因 
为 秦岭 山地 地 处 内 陆 ,冬季 降水 量 少 , 降 水 以 雪 的 
形式 出 现 , 导 致 地 表 反 射 率 增 大 ,融雪 消耗 热量 使 
得 地 表 温 度 降 低 2 ,从 而 使 得 植被 春季 物候 期 推 
US ;秦岭 南北 两 侧 秋季 物候 与 8 一 11 月 的 气温 和 降 
水 呈正 相关 关系 ,光照 充足 ,有 利于 植物 进行 光合 
作用 ,而 较 多 的 降水 会 增加 土壤 水 分 含量 ,使 得 植 
被 生长 可 利用 水 资源 增加 ,从 而 延长 植被 生长 季 长 
度 ,使 得 物候 结束 期 推迟 。 

同时 将 研究 区 内 西安 站 地 面 观测 的 物候 数据 
同 遥 感 提取 的 植被 物候 期 进行 了 验证 。 通 过 收集 
西安 物候 观测 站 2000—2018 EE HW KERL .七 叶 
树 和 灯台 树 4 种 落叶 乔木 的 展 叶 盛 期 和 叶 全 变 色 期 
的 物候 数据 ,发 现 4 种 树木 平均 展 叶 盛 期 为 第 89 d, 
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叶 全 变 色 期 为 第 316 d, 与 遥感 提取 的 物候 始 期 和 物 
候 未 期 集中 天 数 基本 相同 ,表明 MODIS13Q1 物候 产 
品 能 够 反映 秦岭 山地 的 森林 物候 特征 ,可 用 于 秦岭 
山地 的 植被 动态 监测 。 由 于 遥感 监测 物候 反映 的 
是 大 尺度 的 整体 植被 物候 变化 期 ,而 地 面 观测 数据 
多 为 单个 物种 或 单 站 点 数据 ,遥感 物候 与 野外 实地 
观测 的 植物 物候 存在 空间 尺度 扩展 问题 ,但 两 者 在 
物候 的 发 生 时 间 仍 有 很 好 的 对 应 关系 。 因 此 ,寻求 
一 种 合理 的 尺度 转换 方法 进行 遥感 物候 数据 和 地 
面 物 候 观 测 数据 的 精度 验证 ,是 下 一 步 研究 的 重点 
方向 。 


4 结论 


“l 


本 文 利用 DL 模型 ,结合 NDVI 数 据 和 地 面 气 
温 、 降 水 资料 ,提取 了 2000 一 2018 年 秦岭 山地 植被 
物候 生长 季 始 期 .末期 和 生长 期 长 度 , 分 析 秦 岭 山 
地 植被 物候 变化 趋势 和 规律 ,探讨 物候 时 空 变化 及 
其 对 气象 因子 的 响应 ,得 出 以 下 结论 : 

(1) 秦岭 山地 物候 始 期 集中 于 第 60~102 d(3 月 
至 4 月 中 旬 ) ,呈现 从 南北 两 侧 低 海拔 地 区 向 中 高 海 
拔 区 域 推 迟 的 规律 ;生长 季 末期 主要 集中 在 第 315~ 
345 d, 生 长 季 长 度 集中 在 第 225~255 d, 基 本 呈现 出 
从 南北 两 侧 至 山脊 线 逐 渐 缩 短 的 规律 。 总 体 来 看 ， 
物候 期 分 布 格局 与 海拔 高 度 和 地 形 保持 了 较 好 的 
一 致 性 ,海拔 越 高 ,物候 始 期 越 晚 ,未 期 越 早 ,生长 

(2) 从 植被 物候 的 年 际 变化 来 看 ,物候 始 期 
62.25% 的 区 域 呈 现 提 前 趋势 ,53.42% 的 区 域 物候 末 
期 主要 呈现 推迟 的 趋势 ,其 中 59.18% 的 区 域 植 被 生 
长 期 长 度 延 长 , 仅 有 5.12% 位 于 中 东部 中 高 海拔 区 
的 部 分 区 域 显著 延长 。 

(3) 植被 物候 始 期 大 部 分 区 域 与 2 一 5 月 平均 
气温 和 降水 呈 负 相关 关系 ,分 别 占 总 面积 的 50.29% 
Fil 65.24% ; 物候 未 期 与 8 一 11 月 平均 气温 和 月 降水 
量 主要 呈正 相关 关系 ,分 别 占 总 区 域 的 66.63% 和 
59.77%. 

(4) ERITA RRR AERE 4 
季 和 前 期 气象 因子 的 共同 影响 ,不 同 坡 向 物候 变化 
对 气象 因子 的 响应 程度 有 所 差异 ,其 中 秦岭 北 坡 物 
候 始 期 主要 受 当 年 4 月 气温 和 上 年 12 月 降水 影响 ， 
而 秦岭 南 坡 主要 受 当年 3 一 4 月 气温 和 上 年 11 月 降 


水 影响 ; 北 坡 植被 生长 结束 时 间 主 要 受 当 年 9 一 10 
月 气温 和 4 月 降水 影响 ,而 南 坡 主要 受 当 年 9 一 10 
月 气温 和 5 月 降水 影响 。 
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Phenology of forest vegetation and its response to climate 
change in the Qinling Mountains 


YANG Qi, LI Shuheng™, LI Jiahao'’, WANG Jiachuan'” 


(1. College of Urban and Environment Science, Northwest University, Xi "an 710127, Shaanxi, China; 2. Shaanxi Key 
Laboratory of Earth Surface System and Environmental Carrying Capacity, Xi an 710127, Shaanxi, China) 


Abstract: It is important to investigate the vegetation phenology pattern and the driving mechanism of vegetation 
phenology change in a mountain transition zone. In this study, we used the double logistic curve fitting method to 
extract the phenological parameters of forest vegetation in the Qinling Mountains on the basis of the moderate 
resolution imaging spectro radiometer and normalized difference vegetation index time-series images from 2000 
to 2018. Combining temperature and precipitation data, we performed trend analysis, correlation analysis, and 
redundancy analysis to study the characteristics of time and space changes in different phenological periods and 
the response to climate elements in the Qinling area. The following results were obtained: (1) The start and end of 
the growing season ranged from 60 days to 102 days and from 315 days to 345 days, respectively. The length of 
the growth season ranged from 225 days to 255 days. From a spatial perspective, it has evident vertical zone 
distribution characteristics. With an increase in altitude, the start, end, and growth periods of phenology were 
delayed, advanced, and shortened, respectively. (2) The beginning of phenology was advanced, whereas the end 
of phenology was postponed. The proportions of pixels in advance and delayed phenologies were 62.25% and 
53.42%, respectively. The length of the growing season exhibited a lengthened trend of 59.18%. Of the significantly 
extended area, 5.12% area was mainly located in the middle and high altitude areas of the central and eastern 
regions. (3) The initial period of vegetation phenology was negatively correlated with the average temperature 
and precipitation from February to May, accounting for 50.29% and 65.24% of the total area, respectively. The 
end of phenology was positively correlated with the average temperature and precipitation from August to 
November, accounting for 66.63% and 59.77% of the total area, respectively. (4) Redundancy analysis results 
show that the vegetation phenology in the Qinling Mountains was affected by both the current season and the 
earlier period of meteorological factors. The spring phenology was more affected in the phenology occurrence 
period than in the winter of previous year. Compared with autumn weather factors, temperature and precipitation 
in the early period of phenology exhibit a more significant relationship with the end of phenology. The 
phenological changes of different slopes exhibit different responses to meteorological factors. 


Keywords: phenological change; meteorological factors; redundancy analys; Qinling 


